Double-Kasner Spacetime: Peculiar Velocities and Cosmic Jets 
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In dynamic spacetimes in which asymmetric gravitational collapse/expansion is 
taking place, the timelike geodesic equation appears to exhibit an interesting prop- 
erty: Relative to the collapsing configuration, free test particles undergo gravitational 
"acceleration" and form a double-jet configuration parallel to the axis of collapse. 
We illustrate this aspect of peculiar motion in simple spatially homogeneous cos- 
mological models such as the Kasner spacetime. To estimate the effect of spatial 
inhomogeneities on cosmic jets, timelike geodesies in the Ricci-flat double-Kasner 
spacetime are studied in detail. While spatial inhomogeneities can significantly 
modify the structure of cosmic jets, we find that under favorable conditions the 
double-jet pattern can initially persist over a finite period of time for sufficiently 
small inhomogeneities. 
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I. INTRODUCTION 

Astrophysical jets are generally associated with cosmic systems that have undergone 
some form of gravitational collapse. In fact, the direction of the double-jet configuration 
of high-energy outflows is commonly presumed to coincide with the axis of rotation of the 
central gravitationally collapsed source. Similar double-jet patterns have been found in the 
recent analytic studies of the motion of free test particles relative to comoving observers 
in certain time-dependent solutions of Einstein's gravitational field equations |2J. In the 
spatially homogeneous Kasner spacetime, for instance, gravitational collapse occurs along 
one spatial axis, while there is expansion along the other two spatial axes. It can be shown 
that free test particles exhibit peculiar acceleration up and down parallel to the collapse 
axis leading to a double-jet configuration. Indeed, in time, peculiar velocities diminish to 
zero along expanding directions and approach the velocity of light along the contracting 
direction |2| . It is therefore worthwhile to investigate the robustness of these results in the 
presence of spatial inhomogeneities as in the double-Kasner spacetime. We find that for 
small inhomogeneities the double-jet feature involving peculiar acceleration indeed persists 
over a finite time interval. Furthermore, in view of the current interest in "dark flow" , this 
work contributes to the theoretical study of bulk flows associated with peculiar motions in 
spatially inhomogeneous and anisotropic cosmological models. 

Consider a spacetime in which the proper distance along one spatial axis decreases to zero 
asymptotically (i.e., for t — > oo), while the proper distance along each of the other two spatial 
axes tends asymptotically to infinity. Let us introduce a family of fundamental observers in 
such a spacetime with four- velocity vector field = (— gtt)~ l ^ 2 ^ L Q and orthonormal tetrad 
field A M ( Q ,) with \^{q) = We consider only positive square roots throughout. It is natural 
to express the four- velocity vector of a free test particle, = dx^/dr, with respect to the 
fundamental observers via u M = u^X^^y Here r is the proper time along the test particle's 
world line. Studying as a function of time, we have found that there are timelike 
geodesies moving up and down parallel to the axis of collapse such that the Lorentz 
factor of the free test particle relative to the observer family, tends to infinity as t — > oo. 
The remarkable gravitational "acceleration" to the speed of light is an observer-independent 
feature of the double-jet configuration. This feature was first demonstrated in a Ricci-flat 
rotating cylindrically symmetric spacetime 



The timelike geodesic equation in the spacetime of Ref. [l| has the structure of gravito- 
magnetic jets; that is, there are special timelike geodesies for a fixed value of the cylindrical 
radial coordinate that propagate along helical paths up and down parallel to the axis of 
cylinder, which is also the axis of rotation for this gravitational field. The speeds of free 
test particles that follow the special geodesies with respect to the fundamental observers 
(that are spatially at rest) tend to the speed of light as t — > oo. Gravitomagnetic jets are 
attractors. 

We emphasize that in Ref. the physical spacetime domain in which gravitomagnetic 
jets occur is free of matter and singularities. The boundaries of this domain are excluded 
due to the inadmissibility of circular cylindrical coordinates. It is most likely that without 
cylindrical symmetry, an analytic treatment of gravitomagnetic jets would not have been 
possible. Indeed, the main result — namely, the existence of a gravitational mechanism for 
the acceleration of free test particles to the speed of light — has nothing to do with the 
physical limitations of cylindrical symmetry. This follows from our recent study of cosmic 
jets in Kasner spacetimes {2]; that is, we show that similar results hold for the standard 
Kasner models of anisotropic cosmology. 

The purely gravitational mechanism for jet formation that we describe in Ref. Q works 
as gravitational collapse is taking place; once a relatively stationary situation is established 
after the collapse, our mechanism ceases to function. We expect that other, basically electro- 
magnetic, mechanisms would then take over to maintain the jet that has just been formed. 

More precisely, the various rays of high-energy particles that our gravitational mechanism 
originally produces may get confined and collimated where magnetic fields are strong. The 
initial double-jet configuration is then sustained over time by the various MHD mechanisms 
that have been discussed in connection with astrophysical jets. The end result very much 
depends on the MHD aspects of the astrophysical environment. Either our initial streams 
of high-energy particles simply disperse via collisions every which way as cosmic rays or are 
confined, collimated and sustained by the MHD environment, or a combination of both. The 
central "engine", in this approach, is whatever mechanism that sustains the jet in a given 
environment. 

In previous work on tidal dynamics js-9], the relative behavior of a congruence of 
geodesies has been studied. Moreover, tidal gravitational acceleration of ultrarelativistic 
particles has been considered in a general context but via geodesic deviation within the 



small-deviation approximation {4- 9]. To study tidal acceleration/deceleration, imagine a 
congruence of timelike geodesies in a gravitational field; then, choosing one path in the con- 
gruence as representing the world line of our reference observer, we can express the motion of 
the other members of the congruence with respect to the local frame of the fiducial observer. 
A naturally invariant description of relative motion can be obtained by establishing a Fermi 
coordinate system in the neighborhood of the reference observer. The equation of motion of 
a nearby geodesic with arbitrary relative velocity (less than c) with respect to the fiducial 
observer is the generalized Jacobi equation. Indeed, the generalized Jacobi equation in a 
Fermi system has been employed to illustrate the concept of "critical speed" given by c/y/s, 
where s = 2 for linearized gravity in Fermi normal coordinates — a review of this approach 
is contained in Ref. 9]- This circumstance is compared in Ref. js] with electrodynamics 
(s = 1), and, as expected, it is found that the corresponding critical speed is c. When the 
relative speed of the neighboring geodesies exceeds the critical speed cj\f2 ~ 0.7c, gravi- 
tational tidal forces can behave in a manner that is contrary to Newtonian expectations. 
This can lead, in certain circumstances, to the tidal acceleration of ultrarelativistic (i.e., 
v > c/v2) particles. The resulting gravitational acceleration mechanism has been exten- 
sively discussed in the stationary exterior field of Kerr black holes in connection with the 
speed of astrophysical jets 4j-[9|. Instead of a single fiducial observer involved in the spatial 
neighborhood of free test particles, however, one can consider a whole family of observers 
associated with the ambient medium that perform pointwise measurements along the path 
of the free particles. That is, instead of tidal acceleration, we are interested here in peculiar 
acceleration, namely, the acceleration of free test particles relative to the comoving observers 
of the ambient medium. Hence, we adopt a completely different approach here; in the present 
treatment, we deal with exact pointwise calculations within the context of certain specific 
dynamic spacetimes. The critical speed turns out to be c in this exact treatment, which 
proves to be formally similar to electrodynamics. 

To see how this comes about, imagine an arbitrary timelike geodesic world line W with 
unit tangent vector uF" = dx^/dr, where r is the proper time along W. We use units such 
that c = 1 in what follows; moreover, the signature of the metric is +2. Let A At ( Q ) be a 
smooth orthonormal tetrad frame field associated with an arbitrary congruence of fiducial 
observers. In particular, along W we can write u M = vS a > A%). Clearly, = u^A^ 
is coordinate independent, but it does depend on the background congruence of fiducial 
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observers. In particular, any "acceleration" or "deceleration" that could be deduced from 
the study of would depend on the background observers, since a geodesic particle has 
indeed a vanishing acceleration vector; that is, = Du^/dr = 0. The covariant derivative 
of along W vanishes; hence, 

The rate of variation of the tetrad frame along W can be expressed as 



(a) 



*(a) W W(fi), (2) 



dr 

where &( a )(p) is antisymmetric due to the orthonormality of the frame field. It follows from 
Eqs. (|XJ) and © that 

This result is formally analogous to the Lorentz force law. Just as electromagnetic acceler- 
ation of charged particles can be investigated using the Lorentz force law, the gravitational 
acceleration of free test particles relative to the congruence of reference observers can in 
principle be studied using Eq. ((3]). At any given event along the geodesic, the other observer 
families that one could imagine are related to the fiducial family by the elements of the 
Lorentz group. It is therefore essential in the study of this type of gravitational "accelera- 
tion" relative to an observer family to focus on observer- independent features. In fact, the 
divergence of u^°\ the Lorentz factor of the free test particle relative to the congruence un- 
der consideration, along the trajectory is the observer-independent feature that is the main 
focus of our work. That is, the measured speed of the test particle approaches the speed of 
light and this is an invariant property that is independent of any observer as a consequence 
of local Lorentz invariance. This divergence occurs in Ref. [l| and is the central feature of 
gravitomagnetic jets. 

In conformity with the standard practice in relativistic cosmology, in this work we consider 
explicitly the motion of free test particles relative to the background comoving observers. 
Moreover, in the simple examples that we study in this paper, instead of solving Eq. (J5J), 
it is more convenient first to solve the geodesic equation for u M and then investigate the 
behavior of 

u («) = u >*\^ :=7(l,v), (4) 



6 



where v is the peculiar velocity of the free test particle relative to the family of fundamental 
observers and 7 is the corresponding Lorentz factor. 



II. PECULIAR VELOCITIES 

To illustrate the behavior of peculiar velocities in spatially homogeneous spacetimes, 
consider a Bianchi type I model with a metric of the form 



ds< 



df + XV + Y 2 dy 2 + Z 2 dz 2 , 



(5) 



where X, Y and Z are functions of time t. This spacetime is spatially homogeneous with 
three spacelike commuting Killing vector fields d x , d y and d z IQj, 111 ]. 
The gravitational field equations are 



P„ 



(6) 



2 

where T^ u is due to the presence of a perfect fluid with density p and pressure P, 

T lu , = {p + P)U li U v + Pg IBt , (7) 

and the cosmological constant is assumed to be zero. It follows that in comoving coordinates 
with = 5^o, we have 
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Detailed discussions o 



the solutions of these equations for dust can be found, for instance, 



in section 5.4 of Ref. 12] and section 12.15 of Ref. 13] . The timelike geodesies, each with 
four- velocity w M = dx^/dr, are such that the components of along Killing vector fields 
are constants; hence, 
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where C\, C2 and C3 are constants of the motion. Moreover, since w M is a unit vector, we 
find, 

where r is assumed to increase with t along the world line. 

When C\ = C2 = C3 = 0, we have = U^, so that these timelike geodesies coincide 
with the fundamental comoving observers. The orthonormal tetrad frame associated with 
particles of the ambient perfect fluid ("fundamental observers") is then given by A M ( ) = 
and 

A" (1) = (0,^,0,0), A" {2) = (0,0,i,0), A^ {3) = (0,0,0,|). (14) 

We are interested in the peculiar velocity of a free test particle with (Ci, C2, C3) 7^ relative 
to the ambient medium. It then follows from Eq. fll]) that 7 = dt/dr, which is given by 
Eq. (H3J, and 

Ci c 2 c 3 

1Vx = X' 1Vy = Y ' 7 ^ 2 = Z' ^ - 1 

Let m > be the mass of a free test particle and = mu^ be its peculiar four- 
momentum. It is a consequence of Eq. (fT5l) that in general along the x axis, Vt x ) oc X -1 , 
and similarly along the y and z axes. As is well known, a similar relation holds in the 
Friedmann-Lemaitre- Robertson- Walker ("FLRW") models and expresses the law of decay 
of peculiar velocities as the universe expands (see Appendix A). Suppose, however, that 
along one axis — say, the x axis — the universe contracts such that X — > 0; then, Vm — > 00. 
It follows that in this case the peculiar speed of the free test particle approaches the speed 
of light. This issue will be treated in this section in a couple of special cases of interest. 
Let us first consider the behavior of timelike geodesies in the standard Kasner metric 

ds * = _^ + ( l\ 2pi dx 2 + (l) 2P2 d y 2 + (-) 2P3 dz 2 , (16) 

P1+P2+P3 =pI+pI+pI = 1. (17) 

This empty universe model is a solution of Eqs. (lS])- (|lip with p = P = 0. Here to is the 
cosmic time at the present epoch. Henceforth, we will measure time in units of to; therefore, 
unless specified otherwise, we will formally set to = 1 in what follows. We assume that 
Pi < P2 <Ps; that is, 

-^<Pi<0, 0<P2<1, 1<P3<1- (18) 



s 



In this standard expanding anisotropic cosmological model [10, 15, 16], Cartesian coordi- 
nates are admissible for t E (0, oo), t — at the cosmological singularity and (— g) 1 ^ 2 = t. 
Integrating dr = dt/'-f along a timelike geodesic in the Kasner spacetime and choosing the 
integration constant such that r = at t = 0, we find that in general as t goes from to 
oo, r monotonically increases from to oo as well. Furthermore, as cosmic time t increases, 
the Kasner universe contracts along the x axis while expanding along the y and z axes. It 
follows that as t — > oo, v x —> C\/\Cx\, v y — > and v z —> 0. However, for free test particles 
with C\ = 0,v— 7-0 as t— )• oo . Thus in general — that is, except for a set of measure 
zero — all timelike geodesies asymptotically form a double-jet configuration, relative to co- 
moving observers, that is parallel to the axis of collapse. Moreover for t — > 0, v x — > and 
7 — > oo. Thus in general peculiar velocities diminish to zero along expanding directions and 
increase up to the speed of light along contracting directions. These results are all the more 
remarkable because — except for the sign of C\ — they do not depend significantly on the ini- 
tial conditions for geodesic motion. In particular, in expanding directions, the magnitudes 
of the initial velocities of the free test particles do not affect the end result 

It is interesting to consider, in Kasner spacetime, null geodesies, each with tangent vector 
= dx^/d^, where ( is an afline parameter along the path. Then, as before, we have that 

t 2pi ^ = iVi, t 2 ^ = N 2 , t 2 ^=N 3 , (19) 
dQ d( d( 

where N\, N 2 and ^3 are constants of the null geodesic motion. From k^k^ = 0, we find 

^ = (jv?r 2 * + N 2 2 r 2 ^ + Nlr 2 ^) 1 ' 2 , (20) 

where t is assumed to increase with ( along the world line. It follows from an explicit 
comparison of Eqs. ( fT9l -(|20l with the corresponding relations for timelike geodesies that 
timelike geodesies approach null geodesies for t — > oo as well as for t — » 0; indeed, these 
features have been discussed in detail in Ref. {2]] in terms of the notions of speed-of-light 
attractor and repellor, respectively. 

Let us next turn to an examination of timelike geodesies in the Einstein-de Sitter model, 
which is a solution of Eqs. (l8"])- (llip for P = and p = l/(6irGt 2 ). The metric is 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (21) 

where a(t) = (t/t ) 2 ^ 3 and to denotes the present cosmic epoch. The universe in this model 
expands from a singular state at t = to t = 00. Similarly, along any timelike geodesic the 
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proper time r goes from r = to r = oo, once the integration constant is so chosen that 
r = at t — 0. The peculiar momentum of a free particle, V(t), pointwise measured by 
a comoving observer that the particle passes at cosmic epoch t, is proportional to l/a(t). 
Thus in time the peculiar velocities monotonically decay, so that the free particles tend to 
a state of rest relative to the Hubble flow for t — > oo, while for t — > 0, a — > and hence 
7 — > oo. 

There is a particularly simple way of illustrating peculiar acceleration in the case of the 



Einstein-de Sitter universe. The idea is to display typica 



particle orbits using the standard 



conformal time diagram representation (see, e.g., Ref. jl7]). The conformal time rj is defined 
by the relation dt = a(t)dr]. Then, t] = St 2 / 3 ^/ 3 , rjo = 3to and the metric can be expressed 
in the form 

ds 2 = a 2 (r])[-dr] 2 + dr 2 + r 2 (d6 2 + sin 2 ## 2 )], (22) 

where a(rj) = (r]/r] Q ) 2 . We can focus on radial motion with respect to a fiducial fundamental 
observer without loss of generality because of isotropy and homogeneity. As follows from 
the above considerations, there is a constant of the motion given by 

C = ^*, (23) 

where C 2 = C\ + C\ + C 2 , so that the special value C = corresponds to the fundamental 
comoving observers. Expressing Eq. ( |T3l) for this case in conformal time and using Eq. ( 1231) 
result in the equation for the orbits in the form 

J 2 = i +6 y, (24, 

where b = C~ lr qQ 2 . This is an elliptic equation which can be reduced to Weierstrass form by 
using the auxiliary variables a = rj 2 and u = br. Then cr(u) = p(u; ti2, h^) is a Weierstrass 
elliptic function which solves the equation 

a'{u) 2 = Aa{uf - h 2 a{u) - h 3 (25) 



with h 2 = —4:b 2 and hs = 0. Representative solutions for 77 (r) = y/a(br) are shown in 

Fig.m 

As can be understood from that figure, the fundamental observer congruence is stable in 
the future direction and it is attracting a space-filling family of orbits. It is also apparent that 
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FIG. 1: Conformal time diagrams of the Einstein-de Sitter model showing free particle orbits. 
Lightlike geodesies would be at 45 degrees. Vertical lines are worldlines of fundamental observers. 
The curved lines are examples of non-comoving free orbits. The behavior as t — > (corresponding to 
T] — > 0) shows the (negative time) acceleration towards the speed of light relative to the fundamental 
observers (or equivalently the homogeneous hypersurfaces). Left panel: Orbits corresponding to 
one value of the parameter b are shown. Right panel: The family of orbits which is attracted to 
the fiducial observer at r = in the limit t — > is shown. These orbits have different values of b 
corresponding to different initial distances from the fiducial observer. 

there is no acceleration in that direction. In the past direction, the situation is essentially 
the opposite, the flow is unstable and there is relative acceleration to the speed of light, 
v — > 1, as discussed above. In particular, in this case, the same family of orbits that is 
attracted in the future direction is repelled in the past direction. 

We note in the examples above that the proper time can be chosen to be zero, say, at t — 0, 
so that the free particle encounters a curvature singularity when the peculiar acceleration 
to the speed of light occurs at a finite proper time along the geodesic. This circumstance 
appears necessary, since otherwise the geodesic could have been extended. The phenomenon 
of peculiar acceleration to the speed of light may possibly be used as a way of characterizing 
certain spacetime singularities. In short, if it occurs at a finite proper time of the relevant 
geodesic, it implies that the geodesic cannot be continued to that limit and beyond. This 
therefore indicates that the spacetime is timelike geodesically incomplete in the terminology 
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of Rcf. 



12| . It is clear that the singularity is nontrivial in the sense that it cannot be removed 



by extending the spacetime in some way. The more detailed nature of such a singularity is 
an interesting question that is beyond the scope of the present investigation. 

We have thus far considered peculiar velocities in simple spatially homogeneous cosmo- 
logical models. In the rest of this paper, we study peculiar velocities in a simple spatially 
inhomogeneous Kasner-like spacetime. It would be interesting to make a general study of 
peculiar velocities in spatially inhomogeneous cosmological models, as cosmological inhomo- 
geneities tend to mimic dark energy — see 



and the references cited therein. 



III. DOUBLE-KASNER SPACETIME 



To investigate the occurrence of cosmic double-jet configurations in the process of gravi- 
tational collapse, it is necessary to consider a realistic scenario involving a collapsing system. 
However, the actual physical situation involving a collapsing configuration of matter is ex- 
tremely complicated to model properly; therefore, it is useful to consider exact solutions 
of Einstein's field equations that exhibit this property. Imagine, for simplicity, the Kasner 
metric ( JTB"]) such that the metric coefficients also involve functions of spatial coordinates 
(x, y, z). While in this case, at every point in space, simultaneous collapse and expansion do 
occur along the spatial axes, they do so in a spatially inhomogeneous manner. This has the 
consequence that along a timelike geodesic in such a spacetime, as t varies monotonically 
with proper time r by assumption, the peculiar motion involves metric functions that con- 
tain x(r), y(r) and z(t) as well. The influence of spatial inhomogeneities may then totally 
dominate the collapse/expansion scenario. We show that this is indeed the case in certain 
situations via a detailed examination of the geodesies of double-Kasner spacetime. 

Consider, for simplicity, the diagonal spacetime metric 

ds z = - e ^dt 2 + g\dx 2 + g 2 2 dy 2 + g 2 3 dz 2 , (26) 

where go is a constant parameter and for i — 1, 2, 3, we define 

9i ■= e * x t Pi , (27) 

so that spatial inhomogeneity is introduced here via exponential dependence of the metric 
coefficients upon the x coordinate. Here Pi are constants satisfying Eq. (fTT|) as before and 
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qi are given by 

9o + 92 + 93 = 9i, (28) 

9o + 92 + 9s = 9?, (29) 

9o(P2 + ps) = 92(P2 - Pi) + 93(P3 - Pi)- (30) 

The coordinates (i, x, y, z) in this double- Kasner spacetime satisfy Lichnerowicz admissibility 
conditions once t 6 (0, oo) and each spatial coordinate takes values in the open interval 
(—00,00). There is a cosmo logical curvature singularity at t = as in Kasner spacetime and 
(— g) 1 / 2 = texp (2qix); moreover, the spacetime contains two spacelike commuting Killing 
vector fields d y and d z . The connection coefficients for metric (1261) are given in Appendix B. 

It is possible to show that this spacetime is Ricci-flat. Following the work of Kasner fl^ . 
it is natural to look for solutions of Einstein's field equations such that the metric tensor 
depends only on two coordinates, such as, for instance, t and z. Assuming separability, such 
two- variable Kasner- type Ricci-flat solutions were found and discussed in Ref. (19I ]. In our 
case, the metric tensor depends on separable functions of t and x and is thus related to 



solutions of Ref. 



191 ]; in fact, it belongs to Case B of Harris - Zund solutions [19J. Moreover, 



Harris and Zund showed that for normal par ameters such solutions are algebraically general 
(Petrov type I), just as in the Kasner case 19J. The present Harris - Zund solution, essent ially 
in the form of equations (T26T) — (T30l) . was later given explicitly in Appendix 1 of Ref. [ 20] . 
Further discussion of such multi-variable solutions is contained in section 17.3 of Ref. [R 

It is a simple consequence of the algebraic relations (T2~8|) — (T3Q|) that if q = 0, then the 
spacetime is either flat or Kasner. Moreover, if q± = 0, then metric ( 1261) simply reduces 
to the standard Kasner metric. Thus we assume in what follows that go 7^ and q\ ^ 0. 
Now let g = 90/91, 92 = 92/91 and g 3 = 93/91; then, it follows from Eqs. ([28]) and fl29|) 
that the constants q := (go, 92, 93) satisfy the same algebraic relations given in Eq. ( Tl7|) for 
p := (pi,P2,Ps) and can be parameterized in the same way. The relationship between p and 
q, as discussed in detail in Appendix C, is marked by reciprocity. Once p is chosen, then 
Eqs. ( I28l) -(l30l) determine q, thus leaving the metric still dependent upon one free parameter 
qi, which could be positive or negative. Now under the parity transformation x — > —x, 
the metric remains invariant if the g's are replaced by — 9's, which leaves q invariant; to 
remove the ambiguity in the sign of q\ with no loss in generality, we assume henceforth that 
qi > 0. Moreover, we note that metric (1261) is invariant under the exchange of y with z if 
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simultaneously we exchange p 2 with p 3 and q 2 with q 3 . This exchange invariance will be 
used in what follows. 

It is useful to introduce a dimensionless radial coordinate £ defined by 

e qix = £, (31) 



so that as x G (—00,00), we have £ G (0,oo). Under such a transformation, metric ( t26|) 
takes the form 

cfe 2 = -£ 2 *dt 2 + i-t 2pi rf£ 2 + i 24 H 2p2 dy 2 + £ 2q H 2ps dz 2 . (32) 

We recall that Kasner metric (fT6|) is the standard timelike form of the Kasner solution, while 
its spacelike form may be written as [l^] 

ds 2 = -x 2pi dt 2 + dx 2 + x 2p2 dy 2 + x 2pz dz 2 . (33) 

Here, Cartesian coordinates are admissible for x G (0, 00), x = is a curvature singu- 
larity and (—g) 1 / 2 = x. It follows that double-Kasner metric (1321) is indeed a nonlinear 
superposition of the timelike and spacelike forms of the Kasner metric. In the form (l3"2l) . 



the double-Kasner metric belongs to Case A of the Harris - Zund classification 
note that certain other "double-Kasner" metrics have also been considered in Ref. 
Ref. 



91 We 
2l| and 



22J ; in the former, in spacetimes of lower symmetry and in the latter, in 5D spacetime. 
In the Kasner case, the three parameters cannot all be equal. If two are equal, then either 
p = (1, 0, 0), in which case the spacetime is flat or p = (—1/3, 2/3, 2/3), which corresponds 
to the non-flat plane-symmetric case. Otherwise, the three parameters are all different and 
take values within the open interval (—1/3, 1), with one negative and two positive. Let us 
now consider in a similar way metric ( 132]) by treating the special cases of p and q. We first 
notice that p and q cannot be equal in curved spacetime. That is, p = q leads to p\ = 1 or 
qo = 1 or indeed both, since with p = (1,0,0), we find q = (1,0,0); then, in terms of the 
radial coordinate £ defined in Eq. ( 13TT) . the double-Kasner metric reduces in this case to 



ds 2 = -i 2 dt 2 + \t 2 de + dy 2 + dz 2 , (34) 
Qi 

which turns out to represent flat spacetime. Let us briefly digress here and mention a 
generalization of this case such that p\ = qo] indeed, it can be shown that for a non-flat 
spacetime, we must in general assume that pi^ qo- For instance, a double-Kasner spacetime 
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with p = (0, 1, 0) and q = (0, 0, 1) is flat. The other five special cases are not flat and can 
be expressed as: 

I. p= (0,1,0), q=(||-i), (35) 
/ 1 2 2 x - ,3 6 2 N 

2- P=(-g,g, g), Q = , ^ , ~ ^ , (36) 

3- P = (3,-3,3), q=(^,-^,^), (37) 

4. p = (-,-,--), q =(--,-,-), (38) 

,6 2 3, A ,2 1 2 N 

5. p = (-,--,-), q =(-,--,-)• (39) 

Here we have taken into account the exchange invariance as well as the fact that go 7^ and 
that for a non-flat spacetime g ^ q±. 

Finally, let us note that there is a homothetic Killing vector field here just as in Kasner 
spacetime. We recall that in the non-flat case, pi ^ 1 and q ^ q±. Consider the coordinate 
transformation (t, x, y, z) — > (t, x, y, z) given by 

t = e^°t, x = x + u u y = e^y, z = e" 3 z, (40) 

where uq is a constant parameter and 

1 — p\ 

coi = w , (41) 

Qi - qo 

^2 = (qo - <?2 Vl + (1 - £>2 Vo, (42) 

w 3 = (go - <?3 Vi + (l - p 3 W (43) 

Then metric ( l26i) changes only by a constant factor, namely, 

ds 2 = e 2{ulo+qouJl) ds 2 . (44) 
The homothetic generator can be easily deduced from the above considerations. 

IV. SINGULARITIES OF THE DOUBLE-KASNER SPACETIME 

A Ricci-flat spacetime has four algebraically independent scalar polynomial curvature 



invariants. They can be expressed as (cf. chapter 9 of 10j) 



Xi = R^poR 1 """ - iRpu P aR*^ po (45) 
T-% = R^ pa R paal3 R a f3 fJjU + iRfj, U paR paal3 R* a^ v ■ (46) 
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For the double-Kasner spacetime represented by metric (1261) . T\ and X 2 are both real and 
are given by 

X 1 = -IQpmp^C^ - lQq oqi q2q 3 t-^C A ~ Sq^^C 2 ^ K(p, q), (47) 

X 2 = 48p2pgpgr 6 r fl * - 48g 2 g 2 2 g 3 2 t- 6 ^r 6 

+ 24g^- 2 ( 1+2 ^)r 2(2+<? " o) ^i(p, q) - 24g 2 r 2 ( 2+Pl )r 2(1+2(fo) £2(p, q). (48) 

Here we have employed the radial variable £ defined in Eq. fl3TT) ; moreover, to simplify 
the above invariants for the double-Kasner spacetime, we have made extensive use of the 
relations given in Appendix C. The expressions for K, L\ and L 2 are in general complicated 
and are given in Appendix D. 

Double-Kasner spacetime is a nonlinear superposition of two different Kasner spacetimes 
that are usually characterized as timelike and spacelike — cf. Eq. fl32|) . where we temporarily 
set q\ — \. We note that for the standard timelike Kasner spacetime, 

X 1 = -16 PlP2 p 3 t-\ X 2 = 48p 2 p 2 p 2 r 6 . (49) 

Similarly, for the spacelike Kasner spacetime, 

Xx = -\Qq q 2 W\ ?2 = -48g 2 g 2 2 g 3 2 r 6 - (50) 

These invariants can be recognized as the components in the "nonlinear superpositions" 
given in Eqs. (jUj) and (jl8"j) . 

It follows from these curvature invariants that the double-Kasner spacetime has curvature 
singularities at t — and £ = 0, as expected. These can be generally characterized as 
spacelike and timelike, respectively. Moreover, it follows from Eqs. ( )47j) and (H81 that for 
Pi < 0, there is an additional curvature singularity at t = oo; similarly, for q < 0, there 
is an additional curvature singularity at £ = oo. It is important to note that the physical 
spacetime domain under consideration in this paper is free of any singularities; in fact, 
curvature singularities occur at the boundaries of the admissible intervals of t and £. 

V. TIMELIKE GEODESICS 

We consider a free test particle with four- velocity vector = dx^/dr as before. The 
component of this vector along a Killing vector field is a constant of the motion; there- 
fore, g\dyjdr = C 2 and g^dz/dr = C 3 , where C 2 and C 3 are integrals of geodesic motion. 
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Moreover, we find from u^u 11 



-1 that 



;,'/()■'' 



dt 
dr 



C3V/2 



(51) 



where we have assumed that t increases with r along the geodesic world line. The geodesic 
equation for x can be written as 

^ + 2 J + (<?i + *)(^) 2 = [ - <?o + (92 - 9o)Sr + (93 - <?o)Sf] • (52) 

Here we have used Eq. (l5Tj) as well as the connection coefficients for metric ([26]) given in 
Appendix B. It is useful to define a new function W given by 

.da; 



_ f2pi e (qi+qo)x_ 



dr 



Then the autonomous system of equations that we must numerically integrate is 

,W 2 



dt 
dr~ 



-q x 



I _|_ e 2 (<?i-9oK 



9l 



CI 
9l 



CflV2 

9i 



dx 
dr 



m-*>YL d JL-9i d A-9i 



1 ' 

91 



dr 



g\ ' dr 



1 > 

93 



dW 
dr 



c 2 c 2 

Qo + (92 - 9o)-f - + (?3 - <?o N 
#2 



Let us now imagine the fundamental observers in this spacetime with tetrad field A' 



where for i — 1,2, 3, 



A' 



(0 



9. 



■Pi, A> } 



(53) 

(54) 

(55) 
(56) 

(a), 

(57) 



It follows that the Lorentz factor of the free particle 7 and its velocity v as measured by the 
fundamental observers are given in terms of the system ([54 l) -(l56 l) by 

dt 



7 



and 



dx 



dr 



dy 



dz 



iv x = gi — , <yv y = g 2 — , ^v z = g 3 —. 
dr dr dr 



(58) 



(59) 



If the q's all vanish, W — C\ is a constant and system (j5"4"|) - (jS"6"|) reduces to the timelike 
geodesic equation for the standard Kasner spacetime, as expected. In that case, one of the 
p's is negative ("collapse") and the other two are positive ("expansion"), resulting in an 
asymptotic double-jet pattern along the direction of collapse. However, the situation is gen- 
erally different for the spatially inhomogeneous double-Kasner spacetime with nonvanishing 
q's as demonstrated by the following numerical results. 
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VI. NUMERICAL RESULTS 

In integrating the system of equations (!54j) - (l56|) . we must first fix the p's, g's and the 
specific momenta along the y and z directions given respectively by C2 and C3. Starting 
from initial conditions at r = specified by t(0) — to — 1, x(0) = xo, y(0) = 0, z(0) = and 
W(0) = Wq, we numerically integrate forward in proper time such that t — > 00 and backward 
to t — 0. Let us note that in our dynamical system, the equations for t(r), x(r) and W(r) are 
the primary coupled ordinary differential equations. This is due to translational invariance 
along the y and z directions; moreover, we can always choose the initial values of y and z 
at the origin of (y, z) plane without any loss in generality. 

There are an infinite number of possible choices of p's and g's; once p is chosen, the two 
possible q vectors can be algebraically determined, as discussed in Appendix C. We then 
define q = q e, q\ = e, q 2 = q\e and q 3 = q 3 e. Here, e is dimensionless, since e = ct <?i and 
c = to = 1 in accordance with our Kasner conventions; in fact, e G [0, 00) is such that for 
e = we recover the standard timelike Kasner metric. Thus e is a measure of how close 
the double-Kasner metric is to the original Kasner metric; moreover, e _1 is a parameter 
that represents the extent of spatial inhomogeneities in the x direction. It is important 
for the theory of cosmic double-jet configurations to recognize that for sufficiently small 
e, < e << 1, general continuity arguments connect the double-Kasner geodesies to the 
Kasner geodesies over finite intervals of proper time. The primary autonomous differential 
equations can be expressed as 

dt 



dr q 



1 + wH- 2pi c 2q '° + cfr 2p2 «r 2(f2 + c*3 2 t- 2p3 r 2(f3 1/2 , (60) 



= eWr 2pl C 4 °, (61) 

dr 

4o + (% - g 2 )C 2 2 r 2p2 r 2(f2 + (4o - <?3)^" 2p3 r 2(f3 l ; (62) 



dr 

furthermore, dy/dr = C 2 i" 2p2 £" 2 ' f2 and dzjdr = C 3 t~ 2p3 ^" 2 * . For W ^ 0, the primary 
differential equations can also be written as 

dt t Pl £- q '°(W 2 +U) 1 ' 2 



d£ eW 



(63) 



± w2 = JWlO (64 ) 
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where U is given by 

U(t, = t 2pi £ 24o {l + C%t- 2p2 C 2q ' 2 + C 3 2 t- 2p3 r 29 " 3 )- (65) 

Equation (1641) can be simplified in certain situations, thereby leading to the application of 
the effective potential method in one- dimensional motion — see, for instance, Eqs. fl67|) and 
(J68l) in Case 1 below. 

The autonomous system of ordinary differential equations for timelike geodesies has sin- 
gularities at t — and £ = 0. The family of solutions of this system depends continuously on 
proper time r as well as the parameter e. These solutions converge uniformly to the solution 
with e = on every finite time interval. That is, once a finite interval of proper time is 
chosen, then for a sufficiently small e, the geodesies of double-Kasner spacetime behave as 
the geodesies of Kasner spacetime at least up to the end of this fixed time interval. 

We are interested in the main characteristics of peculiar motion in the double-Kasner 
spacetime; therefore, it is useful to mention for future reference that in terms of the radial 
coordinate £, 

iv x = Wt^C 40 , iv y = C^r 9 " 2 , iv z = C 3 t-^r 43 - (66) 

A complete numerical analysis of timelike geodesic motion in the double-Kasner spacetime 
is beyond the scope of our investigation. To simplify matters, we therefore limit our attention 
to the five special cases listed in Eqs. fl35|) - fl39|) . which will be treated in turn below. We hope 
that our approach captures the main features of peculiar motion that could be of interest 
in connection with cosmic jets. 

A. Case 1: p= (0,1,0) and q= (§,§,-±) 

The temporal dependence of the metric in this special case simply consists of expansion 
along the y axis; moreover, the corresponding curvature invariants considered in Sec. IV 
essentially reduce to those of the spacelike form of the Kasner metric. In fact, theses curva- 
ture scalars are given in Case 1 by q\Xi and qfT 2 , where Zi and I2 are the spacelike Kasner 
invariants given in Eq. ( 150]) . 

Numerical experiments involving forward integration in time indicate in this case the 
dominant attractive character of the timelike singularity at £ = corresponding to x = —00. 
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In finite proper time, all timelike geodesies appear to approach this singularity with a peculiar 
velocity that approaches the velocity of light. To see this analytically, let us find dW/d£ 
from Eqs. ( 16TT) and ( 1621) . In this case, the result is 

- = |{i + d(. (67) 

Numerical experiments indicate that £ = at t = t s > 1, r = r s > and = W s < 0. 
Integrating Eq. (}6T]) . we find 

W 2 + (d+C3 2 a = ^ 2 - (68) 

It follows that if Wo > at £(0) > 0, then a free test particle will move such that its x 
component of motion monotonically increases and W monotonically decreases until x, or 
equivalently £, reaches its maximum value where W = 0. At this point, the test particle 
stops along the x axis and changes direction such that subsequently W < and the particle 
falls toward the timelike singularity at £ = 0. It then follows from Eq. (1661) that 

7 = rSfa7 + ^V (69) 



and as £ — > 0, 7 — > 00 and 



*JWJ C 2 



(C| + W)^' ^ (C| +^2)1/2' 

so that u — )• 1. 

Using Eq. (1681) . it is possible to reduce the differential equations for the world line to 
quadratures in this case. It then follows from these results in the case of backward integration 
in time that as t — > 0, £ approaches a finite nonzero number and 

C 

V x 0, Vy~¥ T77T, u 2 -> 0. (71) 

I ^2 1 

Indeed, it is clear from Eq. ( )69|) in this case that 7 — > 00 as t — > 0. It is important 
to recall here that the results for geodesic motion in the spatially inhomogeneous double- 
Kasner spacetime are not in general expected to be consistent with the general notion, 
developed for spatially homogeneous spacetimes in Sec. II, that peculiar velocities decrease 
to zero along expanding directions and increase up to the speed of light along contracting 
directions. However, throughout this section, our numerical results for backward integration 
in time generally agree with those in Kasner spacetime, since £ generally approaches a finite 
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nonzero value as t — > 0. To see the influence of spatial inhomogeneities on peculiar velocities, 
we henceforth concentrate on forward temporal integration. 

These analytic results illustrate a feature that is ubiquitous in the numerical results of 
forward integration in the other well, except for Case 4, where go < 0. Therefore, 

let us assume that % > and we set C2 = C3 = for the sake of simplicity. As we integrate 
forward from r = 0, we soon encounter the £ = singularity at t = t s > 1, r = r s > and 
W = W s < 0. Thus near £ = 0, under the square root in Eq. (160]) . we can neglect the first 
term (i.e., unity) in comparison to the second; hence, 

^ « -wr pi c 24 °, (72) 

GST 

^ = eWt- 2pi C 4 °- (73) 
dr 



Thus 



^ « -er pi £* (74) 
Integrating this relation near the singularity results in 

-&)*;*(*.-*)• (75) 

Similarly consistent results can be obtained for r and W near the singularity. From Eq. (1661) . 
we find that as £ — > 0, 7 diverges; indeed, near £ = 0, 

7 « (-W.^jr*, (76) 

so that v x — > —1, since f y = and v z = 0. 

For the sake of completeness, it is important to compare and contrast here the timelike 
motion of a test particle along the x direction, discussed above, with that of a light ray. 
The equation of motion of the light ray follows from Eq. (|32|) and ds 2 = 0, namely, d£/dt = 
±et~ Pl £ q ~°. We find that if the light ray is initially moving along the positive x direction, 
then integrating the equation of motion with the upper sign implies that it will continue to 
do so and x — > 00 as t — > 00. On the other hand, if the light ray is initially moving in the 
negative x direction, then integrating the equation of motion with the lower sign implies 
that it reaches the singularity at x = —00 in a finite interval of time, essentially as in the 
treatment presented above. 
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B. Case 2: p= (-^, f , f ) and q= (f,|,-f) 

This case is particularly interesting as the direction of collapse coincides with the direc- 
tion of spatial inhomogeneities in the double-Kasner spacetime. Thus for sufficiently small 
e, we expect that test particles following timelike geodesies would exhibit a double-jet con- 
figuration parallel to the x axis as in Kasner spacetime. This turns out to be true only if we 
limit our considerations to finite intervals of proper time. Beyond that, the singularity at 
£ = is expected to dominate the motion as discussed in Case 1. This duality has interesting 
consequences to which we now turn. 

Let us first recall that for e = 0, free test particles with nonzero peculiar velocity form an 
asymptotic Kasner double-jet pattern parallel to the x axis as t — > oo. To simplify matters, 
we assume that C-i = C3 = 0, so that we have geodesic motion only in the x direction. 
Moreover, let us suppose that e, < e << 1, is fixed. Then, a free test particle with a 
positive peculiar velocity characterized by Wq > is expected to move along the positive x 
direction and experience peculiar acceleration just as in the Kasner spacetime. Numerical 
experiments show that this is indeed the illustrated in figure 2, but that peculiar 

acceleration along the positive x axis later changes to deceleration. The particle decelerates 
for a while until it stops (7 = 1) at a finite proper time and reverses course, moving along 
the negative x direction until it eventually reaches the singularity x = —00 (£ = 0) after a 
long but finite proper time with its peculiar speed approaching unity (7 — > 00). To see how 
this comes about, we note that Eq. ( 1621) can be written in this 

dw 3 , 4 

IF = ~7 £r7 - (77) 

Let us define a new parameter r' along the geodesic world line such that dr' = ^~ 4j/7 dT, so 
that t' increases with proper time r along the path. It follows from Eq. ( 1771) that —dW/ dr 1 is 
a positive constant that is much less than unity; therefore, W monotonically decreases along 
the path. This means that even when W is initially positive, Wq > 0, it will eventually turn 
negative in finite proper time and then the general argument presented in Eqs. (172|) - (!76|) 
can be employed to show that the free test particle should stop at some time and then fall 
back toward the dominant singularity at x = —00. An example of this general behavior is 
given in Fig. 2 for Wq = 1. If Wq is negative, the corresponding test particle simply moves 
along the negative x direction and inevitably ends up at the x = —00 (£ = 0) singularity 
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with its peculiar velocity approaching the velocity of light, in accordance with the analysis 
contained in Eqs. fl72|) - (l76|) . 

When we allow C2 and/or C3 to be nonzero, the motion is in general more complicated 
due to the extra degrees of freedom; in fact, an interesting oscillatory behavior can in general 
occur that will be discussed in the last part of this section. 

We conclude that even when spatial inhomogeneities are turned on, as in the double- 
Kasner spacetime, the Kasner double-jet pattern can survive over a finite time interval 
under favorable circumstances, but it is then significantly modified by the presence of spatial 
inhomogeneities. 

Finally, backward numerical integration in the case under consideration in this subsection, 
namely, xq = 0, Wo = 1 and C2 = C3 = 0, is consistent with the result that as t — > 0, 7 — > 1 
and v x — > 0, while v y = v z = by assumption. That is, the motion is in this case confined 
to the x direction, which expands as t — >■ and, as expected, we find that v x — > 0. The 
same result can be obtained from a straightforward mathematical analysis of our dynamical 
system near t = 0. The geodesic equations of motion are simplified by considering the 
dominant terms for t — > 0; in fact, this can be done for all the subsequent cases discussed in 
this section. 

C. Case 3: p = (§,-§, |) and q=(f,-f,f) 

In this case, the Kasner limit (for e = 0) involves a double-jet pattern along the y axis. 
Following our general approach, we set e = 10~ 3 and numerically integrate the equations of 
motion forward in time with xo = 0, Wo = 0, C2 = ±1 and C3 = 0. As expected, the result 
of the integration is similar to the Kasner case over a certain initial time interval, but later 
the motion along the negative x direction toward the timelike singularity at x = —00 takes 
over and 7 tends to infinity. That is, as £ — > 0, v 2 x + Vy — > 1, while v z = 0. The situation 
in the descent toward the singularity is essentially analogous to the analysis contained in 
Eqs. ( !72|) -( !76|) for go > 0. The numerical results for 7 are presented in Fig. 3. 

When we numerically integrate the equations of motion backward in time, we find that 
as t — y 0, v x — > 1, v y — > 0, while v z = 0. Thus 7—7-00 and the peculiar velocity of a free test 
particle approaches the velocity of light toward the cosmological singularity at t — 0. This 
conclusion is in agreement with a detailed theoretical analysis of the equations of motion in 
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FIG. 2: Plot of the Lorentz factor 7 versus proper time r for the peculiar motion of a free test 
particle along the x direction in Case 2. The initial conditions at r = are xq = and Wq = 1; 
moreover, e = 10 -3 and C2 = C3 = 0. Thus at r = we have 7 = v2 ~ 1.4. The figure shows 
that the initial peculiar acceleration reaches a peak around r = 50 and then turns to deceleration. 
Further numerical integration shows that the particle eventually stops at around r = 50000, reverses 
course and accelerates toward the £ = singularity with 7—7-00. 
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FIG. 3: Plot of the Lorentz factor 7 versus proper time r for the peculiar motion of a free test 
particle in the (x, y) plane in Case 3. Note that initially 7 = v2 ~ 1.4 at £ = 1, but in time 7—7-00 
as the £ = singularity is approached. The initial conditions at r = are xo = and Wq = 0; 
moreover, e = 10 -3 , C2 = ±1 and C3 = 0. 

Case 3 near t — 0. 
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D. Case4: p = (f,f,-f) andq=(-i,|,|) 

This case is qualitatively different from the other cases as go < 0, which means that an 
analysis similar to that contained in Eqs. (|72p - (|76|) for forward integration does not apply 
here; that is, forward integration in time is not dominated by the singularity at £ = 0. The 
Kasner limit (for e = 0) involves a double-jet pattern along the z axis, which we expect 
can persist here for a finite interval of time when e is sufficiently small. Our numerical 
experiments indicate similar qualitative behavior for e ~ 10 -3 as for e ~ 1. Therefore, we 
set e = 1 and integrate the equations of motion forward in time with x = 0, Wq = 0, 
C*2 = and C3 = ±1. The initial peculiar acceleration later turns to deceleration until 
the peculiar velocity of a free test particle approaches a constant terminal velocity vector 
whose magnitude is always less than the speed of light. The corresponding behavior of 7 is 
presented in Fig. 4. If the initial value of W is negative, W < 0, the free test particle's 
initial movement along the negative x direction comes to a halt after a while, the particle 
reverses course and, as before, approaches £ = 00 with a terminal peculiar speed that is less 
than unity. The existence of the turning point in this case appears to imply that there is 
a barrier blocking the particle's access to the £ = singularity. For large £,£—>• 00, the 
asymptotic behavior of the equations of motion can be worked out in this case, and we find 
that t ~ £ 7 / 3 , r ~ £ 2 and W ~ £ 2/3 . 

Let us recall here that when q > 0, the only timelike singularity is the one at £ = 0. This 
strongly attracts free test particles, which can reach £ = in finite values of proper time. 
When g < 0, however, there is an additional timelike singularity at £ = 00. Indeed, our 
numerical results in this case seem to indicate that the new singularity at £ = 00 is dominant 
for go < 0, while the £ = singularity is inactive. Moreover, the £ = 00 singularity seems to 
be somehow weaker in terms of its attractive character: free test particles approach £ = 00 
with proper times that go to infinity and terminal peculiar speeds that are less than the 
speed of light. 

This case brings out a certain generic behavior of timelike geodesies for go < 0. To 
illustrate this point, let us set C2 = C3 = for the sake of simplicity and assume that as 
£ — > 00, the asymptotic behaviors of t, r and W are given by 

r~£ K2 , W~C 3 - (78) 
A detailed analysis of the autonomous system f l60|) -f l62|) reveals that this assumption is valid 
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FIG. 4: Plot of the Lorentz factor 7 versus proper time r in Case 4. The initial conditions at r = 
are xo = and Wo = 0; moreover, e = 1, C% = and C3 = ±1. Note that 7 is initially V2 1.4 
and approaches a finite terminal value as r — ?■ 00. 

with K\, k 2 and n 3 all positive and given by 

«1 = — , «2=l+j?l«i, « 3 = <?o +Pi«i- (79) 

1 - pa 

Moreover, it is possible to show that the motion along the x direction reaches a terminal 
peculiar speed as £ — > 00 with Lorentz factor 7^, 

*.=(!- sr. (so) 

The result of the integration of the equations of motion backward in time is that as t — > 0, 
v x — > — 1, v y = and — 0. This means that the peculiar velocity of a free test particle 
approaches the velocity of light and 7 — > 00 as t — > 0. This result is consistent with a 
detailed analytic treatment of the equations of motion in this case near t = 0. 

E. Case 5: p = (§,-§,§) and q =(§, §) 

It is interesting to point out that if in this case we switch p and q, we get Case 3. There 
is a similar connection between Cases 2 and 4. Case 1 stands alone, however, since switching 
p and q in this case leads to Kasner spacetime. 

Let us note that the Kasner limit (e = 0) in this case involves a double-jet pattern along 
the y axis, as in Case 3. Therefore, just as in Case 3, we set e = 10~ 3 and integrate the 
equations of motion forward in time with x = 0, W = 0, C 2 = ±1 and C 3 = 0. Though 
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the details are different, the end result of forward integration is qualitatively the same as in 
Case 3; in fact, this is also the case when we integrate the equations of motion backward in 
time toward t — 0. 

F. Concluding Remarks 

We conclude this section with some general observations. Starting from the present 
cosmic epoch t — to — 1 and integrating the geodesic equations backward to t — 0, we have 
found that the behavior of double-Kasner peculiar velocities near the cosmological singularity 
is essentially the same as in Kasner spacetime. Therefore, we turn to forward integration 
(t : 1 — > oo). Here, beyond the initial Kasner double-jet configuration, our limited numerical 
results for late-time peculiar motion in the double-Kasner spacetime depend significantly 
on whether % is positive or negative. For go > 0, there seems to be a strong source of 
gravitational attraction at the timelike curvature singularity £ = 0. We note that for the 
double-Kasner metric, g u = — £ 2 *\ so that gu — > as £ — > for go > 0. In this limit, 
furthermore, the peculiar velocity approaches the velocity of light. On the other hand, for 
(jo < 0, gu — > as £ — > oo; in this case, £ = oo also happens to be a timelike curvature 
singularity according to the results of Sec. IV. Qualitatively, for % < 0, the £ = singularity 
is somehow inactive and the dominant attractive influence is exerted by the singularity at 
£ = oo, which results in uniform peculiar motion at late times (r — > oo and £ — > oo). 

To find an explanation for this type of late-time behavior, we must turn to the other 
half of the double-Kasner geometry. Let us recall here that the double-Kasner solution 
is a certain nonlinear superposition of the standard timelike Kasner solution (TIBl and the 
spacelike Kasner solution given by 

ds 2 = -£ 2q '°dt 2 + d£ 2 + i 2 ^dy 2 + i 2 ^dz 2 . (81) 

This is a static Ricci-flat solution of general relativity with commuting Killing vector fields dt, 
d y and d z . The admissibility conditions restrict the radial coordinate £ such that £ £ (0, oo); 
moreover, there exists a curvature singularity at £ = 0. The motion of a free test particle 
in this spacetime is such that the components of the four- velocity vector of the particle 
along the Killing vectors are constants of the motion; that is, u ■ dt = —Eq, u ■ d y = E2 and 
u ■ d z = E3. Here Eq > is the constant specific energy of the particle, while, as before, 
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E 2 and E 3 are constant specific momenta. The geodesic equation of motion for the radial 
coordinate £ of the test particle then follows from u ■ u — — 1; that is, 

(|) 2 + V«) = -1, (82) 

where V is the effective potential given by 

7?2 p2 T?2 

V«) = -^| + ^ + ||. (83) 

For (fo > 0, the effective potential is such that the motion can be either confined within 
the interval (0,£ max ] and may then be described as a "fall" toward £ = 0, or it can be 
oscillatory within the interval [£ m in, £ max ]- Here £ min and £ max are turning points such that 
< £min < £max, or the turning points could coincide, in which case £ would be fixed at the 
minimum of the effective potential. Only the fall toward the curvature singularity £ = 
is available for E 2 = E 3 = 0. On the other hand, for go < 0, the effective potential for 
£ : — > oo is monotonically decreasing from +oo to — oo, so that radial motion in £ is 
confined within the interval [£ m i n , oo) and may be described as "escape" to infinity. These 
possibilities for radial motion may be compared and contrasted with geodesic motion in 
exterior Schwarzschild spacetime. Furthermore, we note that d t ■ d t = — £ 2<? °, so that the 
timelike Killing vector becomes null at the endpoints that we have discussed here: the 
curvature singularity £ = for q > and £ = oo for q < 0. 

It remains to discuss the possibility of oscillatory motion within the interval [£ m m, £max] ■ 
Let us assume, for instance, that q\ > % > and g 3 < 0; then, oscillatory motion is possible 
for E 2 0. Moreover, we expect time-dependent oscillatory peculiar motion in the more 
general context of the double-Kasner spacetime. To see an example of this behavior, we 
return to Case 2 above and assume that e = 1, Xq — — 1 and Wo = 0; moreover, we set 
C*2 = 1 and C3 = 2. The time-dependent oscillatory character of the Lorentz factor for 
peculiar motion in this case is illustrated in Fig. 5. We intuitively expect the turning points 
to be time-dependent in this case; in fact, in time £ m ; n approaches the singularity at £ = 
in the case depicted in Fig. 5. 



VII. DISCUSSION 



From the spatially homogeneous examples discussed in Sec. II of this paper, one may 
draw some tentative conclusions (compare, in particular, Fig. [p. As a general remark, 
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FIG. 5: Plot of the oscillatory Lorentz factor 7 versus proper time r in Case 2. The initial 
conditions at r = are xq = — 1 and Wo = 0; moreover, e = 1, C2 = 1 and C3 = 2. Note that 
the Lorentz factor is initially 7 ~ 2.97. The components of the peculiar velocity are all oscillatory 
with increasing amplitudes for v x and v y and decreasing amplitude for v z . In fact, v z tends to zero 
as £ m ; n approaches zero. 

the occurrence of jets requires peculiar acceleration; therefore, jets occur in more special 
situations compared to peculiar acceleration itself, which is consequently a more general 
phenomenon. Restricting attention first to peculiar acceleration, we need some basic ingre- 
dients to describe it, and to characterize its properties. To begin with, we must define an 
appropriate background field of observers. This is because to define peculiar acceleration as 
such, we need some reference background, since only relative accelerations are possible when 
gravitational interactions are involved. In typical physical situations, there is, in most cases, 
a preferred observer family. For example, in spatially homogeneous (non-tilted) models, such 
as Kasner spacetime, the preferred family corresponds to the normals of the homogeneous 
hypersurfaces (see Sec. II). 

As for the observed objects, it is natural to consider them as members of a test field 
of geodesies. It is the peculiar velocity and acceleration of this test field with respect to 
the background observers that is the subject of investigation here. The relative velocity 
can be identified with the peculiar velocity in the astrophysical sense. Therefore, let v 
be the peculiar velocity of a free test particle relative to the background observer family. 
Then our findings indicate that along any expanding axis v — > as expansion tends to 
infinity, while along any contracting axis v —> ±c as contraction tends to zero. Relating this 
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behavior to the background flow, it appears that if in an expanding direction the flow is 
stable (corresponding to attraction) then there is no acceleration. Whereas, in a contracting 
direction, if the flow is unstable (repulsion), then there is acceleration. 

In the gravitational collapse of an astrophysical object such as a star, the peculiar mo- 
tion of free test particles would be referred to the collapsing medium, so that at first sight 
peculiar acceleration might appear as a simple kinematic effect of relative motion. However, 
the limiting situation of peculiar acceleration to the speed of light — emphasized in [1] and 
[2] — is observer-independent and seems to us to indicate true transfer of energy from the 
gravitational field to the free test particles in analogy with the electromagnetic accelera- 
tion of charged particles. Indeed, this gravitational transfer of energy is expected to be a 
general feature of expanding or contracting geodesic congruences in spacetimes that have 
no timelike Killing vector field; in these gravitational fields, energy gain is associated with 
collapse and energy loss is associated with expansion. The energy loss in the context of an 
expanding FLRW model has been discussed by Harrison 23J. We have emphasized energy 
gain in connection with gravitational collapse. However, in realistic collapse scenarios of 
astrophysical interest involving jets and cosmic rays, it remains to see if the gravitational 
energy gain associated with peculiar acceleration is in fact significant. This important prob- 
lem is beyond the scope of our investigation. Instead, we have concentrated in this work on 
the simpler problem of the behavior of peculiar motions in the spatially inhomogeneous and 
anisotropic double-Kasner spacetime. 

The double-Kasner metric reduces to the standard timelike Kasner metric when the 
inhomogeneity parameter e is zero. Thus when e is sufficiently small, < e << 1, we expect 
on the basis of continuity arguments that the Kasner double-jet pattern of peculiar motions 
would persist over a finite time interval. This general behavior is indeed confirmed by our 
numerical work. At late times, however, a significant feature of peculiar motions is that the 
jets tend to move toward either x = —oo (£ = 0) or x = oo (£ = oo) depending on whether 
go > or go < 0, respectively. For go > 0, it follows from forward integration in proper 
time that at late times the dominant aspect of the peculiar motions in the double-Kasner 
spacetime is the strong gravitational attraction toward the timelike curvature singularity 
at x — — oo (£ = 0) with a peculiar velocity approaching the velocity of light. On the 
other hand, for g < 0, we find jets moving toward another timelike curvature singularity at 
x = oo (£ = oo) with uniform terminal speeds that stay well below the speed of light. We 
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have shown that these aspects are associated with the other (spacelike) Kasner geometry 
that is part of the double-Kasner gravitational field. In these considerations the x axis 
is distinguished from the other spatial axes in double-Kasner spacetime, since all spatial 
inhomogeneities occur along the x direction. 

Finally, it is interesting to note that recent measurements of large-scale peculiar velocities 
of clusters of galaxies have been interpreted in terms of an anomalously rapid bulk flow in a 



common direction ( "dark flow" ) — see 



24M27I| and the references cited therein. If confirmed, 



this dark flow would indicate a significant departure from the presumed large-scale spatial 
homogeneity and isotropy of the standard FLRW cosmology. 
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Appendix A: peculiar velocities in FLRW models 

The spatial isotropy and homogeneity of the FLRW universe imply that in the standard 
(t, X-, #5 0) coordinates about any point in space as origin of spherical polar coordinates 
(r, 9, (f>) — with r = RoX, Ro sin \ or R sinh x, respectively, for flat, closed or open models— 
the geodesies issuing from this spatial origin are purely radial. Here Rq is the radius of 
curvature at the present epoch. Therefore, the geodesic equations are such that 9 and are 
fixed while 

— -(l —^ 1/2 
dr \ a 2 



(Al) 



dy D . . 

dr a 2 

where D is a dimensionless constant of integration and a, a(to) = 1, is the scale factor. Thus 
with respect to the natural tetrad frame of the fundamental observer at the spatial origin, 
the nonzero components of Ui a ) are given by um = —dt/dr and uua = D/a. It follows that 

£>V/2 f D 



7 
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Hence as a — > 0, v 2 — > 1. In particular, this proves that in general V(t) oc aft) 1 , which is 



the law of variation of peculiar velocities in standard cosmological models 
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Appendix B: Christoffel Symbols 

The nonzero components of the connection — modulo their symmetry = F^g — for 
metric ( 1261) are given by Y\ x = q , 

rL = jgle- 2 ^, (Bi) 

Ky = f <? 2 V 2 <^, (B2) 
TL = J 9le~ 2qox . (B3) 



Moreover, we have 



1* = * — r^ = ^, (B4) 
9i 1 



* m m = qi, r^ = -g 2 (^) 2 , rL = -g 3 (^) 2 . (B5) 

9i 9i 



Finally, we find that 



r t, = T' T % = ^ ( B6 ) 



and 



r* = j> rL = &- ( B7 ) 



Appendix C: Relations Involving p and q 

We start with the defining properties of p and q, namely, 

J^Pi = l, go + Q2 + g"3 = i, (Ci) 
£p? = i, ?o 2 + g 2 2 + g 3 2 = i. (C2) 

i 

For simplicity, the following consequences of Eqs. fICip and (102[) are stated only for p with 
the understanding that corresponding relations hold for q as well. 
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It follows from squaring Eq. (1C1|) and then using Eq. (1C2j) that 

P1P2 + P1P3 + P2P3 = 0. (C3) 
Next, from Eqs. (ICip and (1C3I) we find that for each i = 1,2, 3, 

R 2 (Pi - 1) =PlP2P3- (C4) 

Squaring Eq. (1C3|) and then using Eq. (IClj) result in 

P1P2 + + = -2piP 2 Ps- (C5) 
Equations ( ICij) and (|C2l) imply that 

J>i =1 + 3 *W 3 . (C6) 

j 

Similarly, it follows from Eq. flC5j) and the square of Eq. ( IC2p that 

^pf = l + 4piP2p 3 . (C7) 

i 

Multiplying Eqs. (1C3|) and f ]C5j) together, and then using Eqs. (1C1|) and (1C6|) . we get 

P1P2 + pIpI + Pbl = ZpIpIpI (C8) 

alternatively, one can cube Eq. (1C3j) . We find from Eqs. (1C3I) . (1C5[) and (1C8I) . via division 

by powers of pxPiPz 7^ 0, that 

E 1 = > E^ = — — . £-4 = —- ( C9 ) 

^ Pi i Pi P1P2P3 i Pi P1P2P3 

Next, multiplying Eqs. (|C2]) and (|C6]) together, and then using Eqs. flC3|) and flCBD . 

we find 

= l + 5PiPaP3- (CIO) 



Similarly, multiplying Eqs. (IC2p and (IU7P together, and then using Eqs. (102[) and flC5p . we 
find 

J^p- = 1 + 6piP 2 p 3 + ?>p\p\v\\ (Cll) 

i 

alternatively, we can square Eq. (IC6[) and then use Eq. (1C8[) . 
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Finally, let us note that the connection between p and q, given by Eq. (l30l) . can be 
written as 

p-q=Pi(l-g ) + go- (C12) 

This relation is invariant under the exchange of p and q, since the right-hand side of 
Eq. (IC12I) can be written as %(1 — p\) + p\. The general reciprocity between p and q is 
noteworthy. Excluding the flat spacetime case with p = q = (1,0,0) given in Eq. (134 j) . let 
us introduce 

P2-P1 n P3-P1 , niQ x 

a = , p = , (013) 

P2 +P3 P2+ P3 

so that Eq. (1GT2I) can be written as 4o = + P43- If is interesting to note that 

(a - f3) 2 + 2 (a + P) = 3. (C14) 

Once one member of the pair (p, q) is chosen, the other can be algebraically determined. 
Suppose, for instance, that p has been fixed; then, one possible q is given by 

go' = ^(2«/3 + a- 3/3 + 3), (C15) 

42 = ^("« 2 + «/? - a + 6), (C16) 
43 = ^(a 2 - ot/3 - 2a + (3 - 1), (C17) 

where A can be expressed as 

A = a/3 -a- /3 + 4. (C18) 
There is, however, a second possible q, which is given by 

4o" = ^( 2 «^ - 3a + /3 + 3), (C19) 
42 = ^-(/3 2 - a/3 + a -2/3-1), (C20) 

gV' = ^(-/3 2 + a/3-P + 6). (C21) 

We note that under the interchange of P2 with p$, or equivalently, of a with /3, the first 
possible q = (q , q 2 , 43) transforms into the second, but with q 2 and q 3 interchanged. 

Let us now turn to a useful Kasner index parameterization due to Lifshitz and 
Khalatnikov — see 
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30] and the references cited therein. It is interesting to extend this 



parameterization to double-Kasner spacetime. Introducing £(u>), 

E(w) = l + w + w 2 , (C22) 
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we define parameters U and V as follows 

1 



U, l + U, U(l + U)\ 



(C23) 
(C24) 



ft=^)(-V, l + V, 7(1 + TO) 

Using this parameterization, Eqs. ( 1C1I) and ( 1C2j) are automatically satisfied. On the other 
hand, Eq. ( lC12p . the relation between p and q, reduces to 



(1 + 2V + UV)(1 + 2U + UV) = 0. 

For a given p, U is fixed and there are two possible solutions for V, namely, 

1 1 + 2U 



(C25) 



V 



2 + U' 



V" 



u 



(C26) 



corresponding to q' and q", respectively. That is, one can directly verify that the substitution 
of V for V in Eq. flC24j) leads to Eqs. (ICT5|) - (ICT7jl . where 



a 



1 + 2U 



U{2 + U) 



(C27) 



Similarly, the substitution of V" for V in Eq. flC24|) leads to Eqs. flCT9|) -( !C2T]) . Explicitly, 



we have 



where 



£ iq '=(2 + [/, (2 + U)(l + U), -l-U), 
S 2 q" =(U(1 + 2U), -U(l + U), (l + f/)(l + 2C/)V 



S x = 3 + 3U + f/ 2 , S 2 = 1 + 3C/ + 3U 2 . 



(C28) 
(C29) 

(C30) 



The same interchange property noted above is recovered in this parameterization when U is 
replaced by 1/U. 

Let us note here, for the sake of concreteness, that with our choice of Kasner p exponents, 
namely, pi < P2 < p 3 , U G [l,oo). In general, there are six possible permutations of the 
Kasner q exponents; for instance, V — > — 1 — V merely interchanges % and q<i in Eq. f ]C24j) . 
while leaving invariant. Such permutations divide the real V axis into six equivalent 
intervals with endpoints that are given by the set {— oo, —2, —1, — |, 0, 1, oo} — see Fig. 4 of 



Rcf. 



30]. 
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Appendix D: K, L\ and L 2 

The purpose of this appendix is to give the expression for K, which appears in the 
Kretschmann scalar X\, as well as the expressions for L\ and L 2 , which appear in the 
expression for I 2 in Sec. IV. 

The curvature tensor can be computed using the double-Kasner metric (126]) and the corre- 
sponding expressions for curvature invariants can be simplified using the relations involving 
p, given explicitly in Appendix C, and similar ones involving q, but not the connection 
between p and q. The results for K, L\ and L 2 are then given by: 

^(p, q) = pI(4o 2 - <M2 + <?2 2 ) + P2(-5o$j + P3<M3) + p 2 (l - go) 

- Pl(l0 2 - PlP2Q2(l + &) - PlP3?3(l + ?s) + Pl?o(l + 2j? 2 ?2 + 2p 3 <? 3 ) 

- PagoSa + pKgo 2 - gbtfa + Q3 2 ), (Dl) 

£i(p, q) = Pi?o(5a 3 + 9s 3 ) - Pife 2 93 2 + P 2 ?o?2 2 (2 - g 2 + 2<f 3 ) + pt%q 3 2 (2 - q 3 + 2g 2 ) 

+ 9o (3pig - 2pi - g ) (P2<?2 2 + P3<?3 2 ) , (D2) 

and 

£a(p, q) = P2P3(3p 2 - 2pi - l)g 2 + p 2 (5pip 2 p 3 + 3piP2 ~ 2pip 3 + 2p 3 )q 2 2 
+ p 3 (5pip 2 p 3 + 3pip 3 - 2pip 2 + 2p 2 )g 3 2 

+ 2p 2 p 3 (l - p 2 g - P 3 <?2 - P 2 gs)- (D3) 

We note that these expressions are exchange invariant, namely, they remain the same under 
the simultaneous exchange of p 2 with p 3 and q 2 with q 3 . 

Next, it is important to implement in Eqs. f ]Dlj) -( lD3|) the relation between p and q given, 
for instance, in Eq. (lC12j) . via the Lifshitz-Khalatnikov parameterization of the double- 
Kasner spacetime described in Appendix C. We recall that for a given p, there are two 
possible values for q, given explicitly by q' and q". The results for the Kretschmann scalar 
can be expressed as 

K(p, q') = -4 Pl p 3 q 'q 2 ', K(p, q") = -4 Pl p 2 q "q 3 " ■ (D4) 
In a similar way, one can show that 

^(p, q') = 2p lP3 q 2 ' 2 q 3 '(3q ' - q 3 ') (D5) 
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and 

L 2 (p, q') = 2p 2 pl(3 Pl - p 2 )%'q 2 '. (D6) 

Moreover, we find that 

£i(p> 4") = 2 PiP2g3 //2 ?2"(3?o" - g 2 ") (D7) 

and 

L 2 (p, q") = 2p 3 pi(3pi - PaKoV- (D8) 

A significant feature of these results is that they are all consistent with the interchange 
property discussed in Appendix C. 
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